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Abstract 

In this article we compute the two-loop SQCD corrections to Higgs-quark-quark couplings in the 
generic MSSM generated by diagrams involving squarks and gluinos. We give analytic results for 
the two-loop contributions in the limit of vanishing external momenta for general SUSY masses 
valid in the MSSM with general flavour-structure. 

Working in the decoupling limit (MsusY ^ v) we resum all chirally enhanced corrections (related 
■ to Higgs-quark-quark couplings) up to order ai"^^^ tan" /?. This resummation allows for a more 

precise determination of the Yukawa coupling and CKM elements of the MSSM superpotential 
necessary for the study of Yukawa coupling unification. 
• The knowledge of the Yukawa couplings of the MSSM superpotential in addition allows us to 

derive the effective Higgs-quark-quark couplings entering FCNC processes. These effective vertices 
can in addition be used for the calculation of Higgs decays into quarks as long as MsusY > 
-^Higgs holds. Furthermore, our calculation is also necessary for consistently including the chirally 
enhanced self-energies contributions into the calculation of FCNC processes in the MSSM beyond 
leading order. 

At two- loop order, we find an enhancement of the SUSY threshold corrections, induced by the 
quark self-energies, of approximately 9% for i^i = MsusY compared to the one-loop result. At 
the same time, the matching scale dependence of the effective Higgs-quark-quark couplings is 
significantly reduced. 
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I. INTRODUCTION 



In the MSSM diagrams with sfermions and gauginos as virtual particles generate impor- 
tant loop corrections to Higgs-quark-quark couplings. After the spontaneous breaking of 
SU{2)l ® U{1)y at the electroweak scale, the Higgs fields acquire their vacuum expectation 
values (vevs) and the genuine vertex corrections to Higgs-quark-quark couplings also gener- 
ate chirality changing quark self-energies (or self- masses) . Thus, there is a one to one corre- 
spondence between loop-corrections to three-point Higgs-quark-quark functions and quark 
self-energies: The correction to a Higgs-quark-quark coupling is given by the corresponding 
chirality-changing self-energy divided by the vev of the involved Higgs field. 

This means that we can simplify the calculation of three-point functions by reducing 
the problem to the calculation of two-point functions (self-energies). In this way, the self- 
energy contributions to quark masses can be directly related to effective Higgs-quark-quark 
couplings which allows for an efficient calculation of the effective Higgs vertices. 

The quark self-energies also modify the relation between the Yukawa couplings of the 
MSSM superpotential and the quark masses (extracted from low-energy observables). Es- 
pecially if tan /3 (the ratio of the vevs of the two Higgs fields) is large, these contributions 
are generically very large and can be of order one [l|,l2|, Isf. In an analogous way, also the 
relation between the CKM matrix of the superpotential and the physical one is altered (by 
chargino-squark diagrams in the MSSM with MFV BS, and in addition by squark-gluino 
diagrams in the general MSSM 0, Isl). Due to these corrections the physical quark masses 
and the measured CKM elements no longer equal the ones which appear in the MSSM su- 
perpotential. One says that theses relations are modified by so called threshold corrections, 
i.e. by the decoupling of heavy particles. Since in Higgs decays, Higgs mediated FCNCs 
(like -Bs(d) mixing and -Bs(d) — ^ At+/i~) and in higgsino vertices the Yukawa couplings (of 
the superpotential) and not the physical quark masses enter, a precise knowledge of these 
quantities and thus of the threshold corrections is necessary. Furthermore, in GUT models 
with Yukawa coupling unification not the effective Yukawa coupling of the SM, but rather 
the Yukawas of the superpotential unify and the SUSY threshold-corrections must be taken 
into account in order to judge if they actually do unify jo], [lo|. In conclusion, it is desirable 
to know the relation between the parameters of the MSSM superpotential and the physical, 
i.e. measurable quantities very precisely. 

Having the relation between the Yukawa couplings (CKM elements) of the superpotential 
and the physical quark masses (physical CKM elements) at hand, one can calculate the 
effective Higgs couplings entering FCNC processes which include the SUSY loop-corrections. 
This is most easily achieved by matching the MSSM on the Two-Higgs-Doublet-Model of 
type three (2HDM HI). The loop- induced couplings of quarks to the "wrong" Higgs field, 
i.e. to the Higgs which is not involved in the Yukawa term in the superpotential, induce 
fiavour-changing neutral Higgs couplings after switching to the physical basis in which the 
quark mass matrices are diagonal in fiavour-space. These effective Higgs couplings can be 
expressed entirely in terms of the physical masses and self-energies depending on MSSM 
parameters. Here a complication arises because these self-energies must be calculated using 
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the Yukawa couplings and the CKM elements of the superpotential which must have been 
determined previously in the process of renormalization by including the loop-corrections, 
i.e. by resumming the threshold corrections. This problem can be solved analytically in the 
decoupling limit of the generic MSSM in which the self-energies are at most linear in the 



Yukawa couplings [U 



The importance of these threshold corrections and thus of the chirally enhanced self- 
energies motivates their calculation at NLO in q;,^. In the MSSM with MFV these corrections 



have been calculated in Refs. 12|, ll3|], 1^ and 15|, ll6|. Here we want to extend this analysis 



to the MSSM with generic sources of flavour violation and resum all chirally enhanced 
effects using the results of Ref. 0, Q, Ell, 17|. In addition, working in the approximation of 
vanishing external momenta, we are able to give relatively simple analytic expressions for 
the self-energies and therefore also the resummation of all chirally enhanced corrections can 
be (and is) performed analytically. 

After discussing the quark self-energies (and their connection to Higgs- quark- quark cou- 
plings in the decoupling limit of the MSSM) in the next section, we derive the relations 
between the MSSM Yukawa couplings and the quark masses at LO in Sec. [ml As the mam 
result of this article we calculate the SQCD contribution to the chirality-changing self-energy 
at the two-loop level in Sec. llVi In Sec. |V]we discuss the topics of Sec. Illll at NLO. In Sec. I VI I 
we derive the effective Higgs-quark-quark couplings and conclude in Sec. IVIII Various ap- 
pendices summarize the relevant one-loop results. 



II. QUARK SELF-ENERGIES, EFFECTIVE LAGRANGIAN AND THE DECOU- 
PLING LIMIT 

As described in the introduction, there is a one to one correspondence between chiral- 
ity changing self-energies and Higgs-quark-quark couplings: in the decoupling-limit of the 
MSSM (MsusY > V and Msusy > Pi where p is the external momentum) chirality changing 
self-energies are proportional to one power of a vev only, and the corrections to the Higgs- 
quark-quark couplings can be obtained by dividing the corresponding self-energy by the vev 
of the Higgs field involved. Thus, as long as the momentum flowing through the Higgs is 
small compared to the SUSY masses and the SUSY masses are heavier than the electroweak 
vev, the decoupling limit is a valid approximation. In this approximation the calculation 
of the Higgs-quark-quark three-point-function can be reduced to the calculation of quark 
self-energies. For this reason we will consider the quark self-energies in this section in some 
detail and discuss the decoupling limit. The analysis is valid independent of the loop-order 
(concerning Ug corrections) at which the self-energies are calculated. 

In general, it is possible to decompose any quark (or any fermion) self-energy into 
chirality-fiipping and chirality- conserving parts in the following way: 

^.iv) = ^(P^) + ^S^f ^(p^)) Pn + {^f\p') + "(P^)) Pl ■ (1) 
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Note that the chirahty-flipping parts xi^^^^'^^ have dimension mass, while the chirahty con- 
serving parts are dimensionless. 

In the following we will be interested in the contributions to Eq. ([T]) which involve heavy 
SUSY particles. The reason for this is that only these contribution lead to the threshold 
corrections entering the relation between the quark masses and the Yukawa couplings of the 
MSSM superpotential. It is convenient to work in an effective field theory in which the part 
of the effective Lagrangian containing mass terms and kinetic terms for the quarks is given 

by 

+ (5,. - cjr) ojr + {5f^ - cj!^') o)^' , (2) 



with the operators defined as 



= QfPLqi , Ojf ^ = iqf0PL Qi , 
^ = qjPRq^ , 0;f ^ = zqJ0Pn • (3) 



Throughout this paper, the Wilson coefficients in the effective Lagrangian (|2]) (or, equiva- 
lently, the operators) are renormalized in the MS scheme. The final results for the Wilson 
coefficients will be written as an expansion in Qs, where Qs is meant to be the MS renor- 
malized strong coupling constant of the effective theory, running with 6 (quark) flavours. 

In Eq. (|2]) the term —VgY^%^ 6fi denotes the part of the Wilson Coefficient of the operator 
Oj,- which is induced at tree-level by the Yukawa coupling of the MSSM superpotential. 
The running of VgY^J^^ (and also that of CJ.^^) is the same as the one of the quark mass in 
the SM (in the MS scheme). At the matching scale m-susY, ^tree is just the Yukawa coupling 
Y'^ of the MSSM superpotential^. Note that Y^^l^ is not the effective Yukawa coupling of the 
SM which instead is obtained from the physical quark mass (see Eq. f fTT]) ). 

The Wilson Coefficients cj^^'^^ and cj^^'^^ in Eq. ([2]) contain the effects of heavy 
particles only. Self-energy diagrams involving no heavy SUSY particles, i.e. ordinary QCD 
corrections containing only quarks and gluons, do not contribute to the Wilson Coefficients 
in the matching procedure, because they are the same on the full side (the MSSM) and on 



^ The matching calculation for i^rcc most easily done by using the MS scheme, both on the MSSM side 
and on the effective theory side. When working up to order ag, wc get at the matching scale ttt-susy: 
Ytrcc — ^"^^ 7 where V^' denotes the Higgs-quark-quark coupling of the MSSM in the MS scheme. However, 
it is well known that one should use the DR-scheme on the MSSM side, such that supcrsymmetry is 
preserved. This can be achieved by the shift = (1 + ^^^^wr- This issue will be considered in more 
detail in Sec. El The matching condition then reads: y^^ee = (1 + ^f)Y^. 
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the effective side (tlie 2HDM III or tlie SM). At tlie matcliing scale msusY we find for tlie 
Wilson Coefficients of Eq. using the results for Sf^^gf'(O) and S^^^^f (0) given in Eq. (|A2|) : 



X? In fx? 



a 



Air 



3-Axl + xi + (4xg - 2xi) In (xg 
2{l-xlf 



> atLOincto 



^ tree 



(4) 



Further, in the following we will focus on the non-decoupling pieces of Eq. ([T]), i.e. those 
contributions which do not vanish in the limit MgusY — ^ C)0 (which also includes /i — )■ 
oo). In contrast, all parts which vanish in this limit are called decoupling. There are two 
different kinds of decoupling contributions concerning self-energies (or effective Higgs-quark 
couplings) : 



The first kind of decoupling effects is related to the expansion of the self-energies in 
powers of p^/Mgugy. This expansion is certainly possible in on-shell configurations 
because the SUSY particles are known to be much heavier than the external quarks. 
In this series, higher order contributions are clearly suppressed for all light quarks and 
even for the top quark non- decoupling corrections are only of the order 'ni^ I ^h^'SY — 
4% with respect to the leading term. Thus, higher orders in p^/M|ugY can be safely 
neglected as long as the external momentum p"^ is small, which is the case for all 
low-energy flavour observables. 

The second kind of decoupling effect is related to the mixing matrices (and also the 
physical masses) of the MSSM particles (squarks and charginos/neutralinos) which 
appear because the mass matrices of the SUSY particles are not diagonal in a weak 
basis. These mixing matrices and mass eigenvalues can be expanded in powers of 
v/MsvsY and also in this case it turns out that the decoupling limit (i.e. the leading 
order w/Mgugv) for realistic values of SUSY masses^ is an excellent approximation to 
the full expressions Beyond the decoupling limit higher dimensional operators 

involving several Higgs fields would appear. 



From dimensional analysis we see that all non-decoupling contributions are contained 

^ — 0. Furthermore, the non-decoupling part of 

0) is linear in v. Thus, in 
J- — 0), Syj^"''"^(p^ = 0) and only keep the 

leading term in v which is equivalent to considering operators up to dimension 4 only. This 



in j^'j^^'^^ and xij^^^'"^^ evaluated at p 
Y,q^RR,LL^^2 _ g-j independent of a vev, while Sl^^''^^(p^ 



the following we will work in the limit T!^rf^'^^{p^ 



^ The new results of the CMS collaboration [iSj and the Atlas experiment 19[ require that squark and 
gluino masses are at least of the order of 1 TeV. 
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simplification allows us to perform a analytic resummation of all chirally enhanced effects 



as developed in Ref. [11 



There is a fundamental difference between Yf-j^^'^^ and Yf^^^'^^ (and thus also between 
QqLR,RL QqRR,LL^ eveu though both pieces do not decouple. We explain this issue at 
one-loop order: Y'jf^'^^ enters always proportional to the quark mass itself into the renor- 
malization of the Yukawa coupling and CKM elements and thus has the same generic size 
as an ordinary QCD loop correction (it is of order a^). Furthermore, as we will see later, 
the even do not contribute to effective Higgs-quark-quark couplings at the one-loop 

level |20| . On the other hand, Yf^^^'^^ can be "chirally enhanced" by a factor tan /3 
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or yl^/(F^^MsusY) 0] which can compensate for the loop-factor. Due to this possible en- 
hancement, Y'jf^'^^ generates the most important contribution to the threshold corrections 
between Yukawa couplings and quark masses. The resulting Wilson Coefficients cj^^^^^ 
can even be of order one, i.e. numerically as large as the corresponding physical quantities 
{rriq. in the flavour conserving case or Vfi x max[mq.,mq.] in the flavour changing one). 
Furthermore, concerning flavour changing neutral Higgs couplings T.'j- ' even constitutes 
the leading order, since these couplings are first generated at the one-loop level. 

Because the gluino contribution to Y'jj^^'^^ involves the strong coupling constant, it is 
the numerically dominant contribution to the threshold corrections modifying the relations 
between the quark-masses and the Yukawa coupling. Regarding flavour changes, in the 
MSSM with MFV only the chargino contribution enters the renormalization of the CKM 
matrix, but once there are sizable non-minimal sources of flavour- violation, again the gluino 
contribution becomes dominant. The neutralino contribution is in most regions of parameter 
space suppressed (except if the gluino is much heavier than the other SUSY particles). Thus 
we consider the gluino contribution in this article. The calculation of the chargino and 
neutralino induced contributions to the threshold corrections and the effective Higgs-quark- 



quark couplings is work in progress 22 



From the arguments given above we see that at any loop-order (concerning corrections) 
the chirality-flipping quark self-energy containing at least one gluino and one squark as 
virtual particles is always proportional to one^ off-diagonal element A?^^^ of the squark 
mass matrix which, in the super-CKM basis, is given by 



(5) 



with A^/"^ = AJP*. Note the presence of the tilde in the Yukawa couplings V^K This refers 
to the fact that a squark-squark-Higgs coupling is involved, while V^^ entering the Wilson 
coefficient Yt^ee Eq. ([2]) is a quark-quark-Higgs coupling. Of course, both of these couplings 
are a priori equal in the MSSM due to supersymmetry, and could be identified with each 



More precisely, in the decoupling limit Sy-.f^ is linear in A'*-'^^, while beyond the decoupling limit it 
contains all add powers of A''-^-'*. 
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other from the beginning if the calculations of the chirality-flipping quark self-energies would 
be performed in the DR-scheme, in which supersymmetry is preserved. However, we decided 
to work out in an intermediate step the SQCD two-loop corrections to the self-energies in 
the MS-scheme, i.e., in dimensional regularization followed by modified minimal subtraction 
rather than using dimensional reduction. At this level, the two couplings and Y^' are 
different and therefore have to be distinguished in the notation. We will discuss this in more 
detail in Sec. |Vl 

The elements A'f^^^ generate chirality-enhanced effects with respect to the tree-level quark 
masses if they involve the large vev Vu (tan /3-enhancement for down-quark) or a trilinear 
AW^-term (4j7(Fi^MsusY)-enhancement). 



A. Decomposition of quark self-energy contributions 

We diagonalize the full 6x6 squark mass matrices in the following way^. 

W^^M^W^ = dmgiml,ml,ml,ml,ml,ml) , (6) 

where rriq^ (s = 1, 6) denote the physical squark masses. 

In the decoupling limit, i.e. to leading order in v/M^vsy, the chirality-flipping elements 
^qLR neglected in the determination of the squark mixing matrices and the 

physical squark masses m?^. The down (up) squark mass matrices are then block-diagonal 
and diagonalized by the mixing matrices T^^^, i^uLy ^ur) following way: 

WHMIWL = diag(m;-.,m^-.,mJ-.,mJ-.,mJ-«,mJ-«), WL = ) ' 

The 3 X 3 - matrices Fg^ and ^qji {Q = U, D) take into account the flavour mixing in the 
left-left and right-right sector of sfermions, respectively. It is further convenient to introduce 
the abbreviations 

kqRR _ -pim -pjm* /q\ 
^^mij — QR'- QR ' \°) 

where i,j,m = 1, 2, 3 and the index m is not summed over. 

On the other hand, left-right-mixing of squarks is not described by a mixing matrix, 
but rather treated perturbatively in the form of two-point qp-qf vertices governed by the 
couplings A^j^^, i.e. by what is called mass insertions 25 . 

For the relations between the Yukawa couplings and the quark masses (to be discussed 
in Sec. IIIip and for the effective Higgs- quark- quark vertices (see Sec. IVl|l it is necessary to 
decompose C^^ ' according to its y^-dependence as 



c: 



dLR 

a 



+ r*. (9) 



Note that our mixing matrices correspond to the hermitean conjugate of the matrices Tq defined in 
Ref. 0,1241. 
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where, as the notation imphes, C'^^^ is independent of a Yukawa couphng. Note that we 
did the decomposition with respect to the Yukawa couphng y*, as Cj^"^ can only involve 
but not (see Eq. (P). 

For the discussion of the effective Higgs-quark-quark vertices in Sec. |VT] we also need a 
decomposition of Sj/"^ and thus of CJ^^^ into its holomorphic and non-holomorphic parts^. 
In the decoupling limit (and in the approximation rrig = 0) all holomorphic self-energies are 
proportional to A-terms. Thus we denote the holomorphic part of the Wilson Coefficient 
as Clf, while the non-holomorphic part (which can be induced by the fi term or by an A' 
term) is denoted as C'^^^^. This means that we have the relation 



III. RELATIONS BETWEEN QUARK MASSES AND YUKAWA COUPLINGS 
AT LEADING ORDER IN 

Let us discuss the renormalization^ of quark masses and Yukawa couplings induced by 
non-decoupling self-energy contributions to the Wilson Coefficients and Cj/'^'^^ 

in the MSSM. For this purpose we focus on the flavour-conserving case but we will return 
to the flavour-changing one in Sec. IVII As it turns out, flavour- changing self-energies only 
contribute to the relation between quark masses and Yukawa couplings at higher orders in 
the perturbative diagonalization of the quark mass matrices. 

For the renormalization and the inclusion of the threshold corrections it is very important 
to distinguish between the Yukawa couplings of the MSSM superpotential and the "effec- 
tive" Yukawa couplings of the SM (or the 2HDM of type III) Y^^ = niq./vq. At the matching 
scale MgusY the running quark mass m^. of the SM is related to the Yukawa coupling of the 
MSSM in the following way: 

^.Kl = m,^ = {v^y^^le + ^) X (l + i (c:^^^ + C^^)) . (11) 

The term - (cf^^^ + C^^^^^ originates from rendering the kinetic terms of the effective 

theory diagonal, or, equivalently in the full theory from the LSZ (Lehmann - Symanzik - 
Zimmermann) factor which originates for the truncation of the external legs. 

As discussed in the last section, only S^j'^'^ (or equivalently Cfj^^ in the effective theory) 
can be chirally enhanced. If we restrict ourselves to this term we recover (in the decou- 
pling limit in which Cfj is proportional to one power of V^' at most) the well-known 
resummation formula for tan /3-enhanced corrections, with an additional correction due to 



^ With (non-)holomorphic we mean that the loop induced Higgs couphng is to the (opposite) same Higgs 
doublet as involved in the corresponding Yukawa coupling of the MSSM superpotential. 

^ Throughout this article, renormalization is not only understood as the process of removing divergences, 
but also as the altering of the relations between different quantities induced by loop contributions. 
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the y4-terms 12|] (and possibly the A'-terms). The resummation formula at leading order is 
given by'' 

^dLR(l) 



Vd 



Vd{l + tan/3e' 



d{i] 



(12) 



with ef'-^'' and C'^.^^^^ defined through Eq. iQ. The superscript (1) denotes the fact that 
corresponding quantity is calculated at the one-loop order. 



IV. CALCULATION OF THE WILSON COEFFICIENT Cjf ^ AT NLO 

In this section we describe the calculation of two- loop contribution to Cj-^^, discuss the 
issue of renormalization, show the expected reduction of the matching scale dependence and 
discuss the decoupling limit in which only one coupling to a vev of a Higgs field is involved. 
In order to be specific, we describe in the following the calculation and the results for down 
quark, i.e. Cj^^, and mention at the very end how C^j^^ can be obtained. 

In the following we write the Wilson Coefficient Cj^"^ as 

QdLR _ QdLR(l) _|_ ^dLR{2) _|_ ^^^-j 

where Cj^^^^^ and Cj^^'"'^^ denote the one- and two-loop contribution, respectively. We 
perform the two-loop matching calculation (order a^) for the Wilson coefficient Cj/"^ in 
D = (4 — 2£:)— dimensions, using dimensional regularization, both for the full theory (MSSM) 
and for the effective theory in Eq. (|2]). The complete list of genuine 1-PI two-loop diagrams 



contributing in the full theory is shown in Fig. [T] (generated with FeynArts [271, |28[). 

As the first two diagrams (involving squark tadpoles) give rise to some subtle points 
concerning renormalization, we ignore them in this subsection and take account their impact 
on Cjj^^ only in the next subsection. 



A. Matching calculation for Cj^ ignoring tadpoles 

In the full theory we first calculate the 1-PI two-loop diagrams (diagram 3 - 16 in Fig. [1]) 
in the approximation niq = and p = 0, but to all orders in v/nisvsY (using exact diago- 
nalization of the squark mass matrices). All diagrams except diagram 16 can be calculated 
by naively setting rrig = and p = 0. Diagram 16, however leads to two contribution: The 
hard contribution, which amounts to the naive limit of vanishing quark masses and external 
momenta of the full two-loop diagram, and the soft contribution which amounts to the same 



^ For large flavour-changing elements also a contribution involving two self-energies can be important for 
the renormalization of the light quark masses [2^. In this case the resummation formula reads for i — 1: 

r'd.LRr'd.LR 



Vd {I + t£in 134) 
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limit but only for the heavy one-loop subdiagram 29|. As the soft contribution is identical 
to the one-loop gluon correction to —i (j'j^^^^'^^ Ojl"^ in the effective theory^, this contribu- 
tion drops out in the matching for C'j^^^^^\ As this soft contribution is the only one which 
is infrared singular, this means in particular that C'j^^^^'^^ is free of infrared problems, as it 
should be. 

We then add the counter-term contributions in the full theory which are induced by 
the renormalization of the parameters m?^, rrig and ag in the corresponding one-loop result 
(where at this level of the calculation these three parameters are renormalized in the MS- 
scheme). The explicit expressions are listed in Sec. lA 3 a[ In one of these counter-term 
contributions the squark-mass counter-term Sm"^^ enters. Of course, when ignoring the 
tadpole diagrams in this section, the tadpole contribution to Sm'^^ also has to be ignored. 

Besides the renormalization of the parameters in the full theory, we also have to attach 
one-loop wave function renormalization constants for the external quark legs to the corre- 
sponding one- loop result. These wave function renormalization constants have two contribu- 
tions: One from a self-energy with a gluon-quark loop and another one from a gluino-squark 
loop. The first one is also present in the effective theory and consequently drops out in the 
determination of C'^^^^'^\ while the second one contributes. Since we perform the renor- 
malization in the MS-scheme only the divergent pieces of S^^^^'^'^ enter C'j^^^'^^ while the 
nmte part gives rise to C^^ 

We now turn to the effective theory. Here, we have to work out one-loop QCD correc- 
tions to -I cj.^^^^''^) Oj,^^, i.e. the IPI dia gram, attach the wave function renormalization 
constants and take into account the effect of the (MS) renormalization constant 5Zo of the 
operator Oj^^. While the first two get canceled against contributions in the full theory (as 
already mentioned above), the effect of the renormaliztion constant of the operator enters 
the matching condition for C'^^^^'^K 

Putting things together, we get the following (schematic) matching equation: 

- 1 6Zo C^^^^^^'-") - t C^',^^^') = Ds + ... + D,, + D\^'' - I [CT„, + CT„-^ + CT„J 



1 

-i- 
2 



^d LR{l,d) _ ^^4^ 



Here CT^g, CT^^^ and CT^^ stand for the contributions induced by the insertions of the 
corresponding counter-terms into the one-loop diagram and Di represents the contribution 
stemming from diagram i of Fig [1] As already mentioned, we did our two-loop calcuation 
in dimensional regularization. So far the parameters rrig, rriq^ and as appearing in the full 
theory were renormalized according to the MS scheme. Also the various Z— factors appearing 
in Eq. (JHj) are renormalized in the MS scheme. The result for C^^'^'^^^^ we get at this level 
corresponds to the sum of the first five terms on the right-hand side of Eq. (1231) . When 



^dLR{i,D) ^j^^ one-loop Wilson coefficient in D = (4 — 2e)— dimensions, i.e. — Y^j^^J^^O), see 

Eq. dH]). 
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giving the explicit expressions for these terms, we freely made use of the unitarity properties 
of the mixing matrices. 

We should be more precise concerning Qs (or as)- In our calculation of the full theory 
side gs stands for gsx, i-e- for the strong couling constant of Yukawa type of the full MSSM 
renormalized in the MS scheme. As we want to express the final result for the Wilson 
Coefficient Cjl'^^'^^ in terms of gflj^, i-e. by the strong coupling constant of the SM in the 



MS scheme running with six flavours, we make use of the relation 30|, |31 



1 + — - 

in 3 



[n, 



6) In(xJ) - J2 (H^s) + HVs)) +^Ca-3Cf 



s=l 



a 



s,MS 



(15) 

Actually, this relation summarizes three steps: First the transition from gg^v in the MS- 
scheme to gs of the full MSSM in the DR-scheme; second the decoupling of the SUSY 
particles, leading to gs running with 6 (quark) flavour in the DR scheme; third the transition 
to gfl^. Eq. (dS]) leads to the additional piece C'j2)>«hiftae 

In principle we should have performed our calculation (of the full theory side) using dimen- 
sional reduction, which preserves supersymmetry, followed by modifled minimal subtraction. 
The corresponding result for Cj^^^'^'^ can be reconstructed by also shifting the parameter 
rriq^ and rrig from the MS— scheme to the DR— scheme in the expression for C'j^^^^^\ As only 
rrig gets such a shift at the relevant order in a^, we denote this contribution in Eq. ( 123|) as 

(2),m5 — >mg 

This completes the derivation of the matching condition for C^j^^*"^'' when ignoring the 
tadpol contribution (i.e. diagrams 1 and 2). Note that we performed our calculation us- 
ing the expression for the gluon-propagator in an arbitrary i?^— gauge and found a gauge- 
invariant result for C'j^^^^'^\ 



B. The squark tadpole 

The diagrams containing a squark-tadpole self-energy as a sub-diagram require close 
examination. Diagram 1 vanishes but the squark-tadpole contained in diagram 2 contains 
a divergence which enforces a renormalization of both the physical squark masses and the 
trilinear couplings of squarks to the Higgs fleld (the Yukawa couplings and the A-terms). 
Thus it has to be decomposed into the corresponding two parts. 

Let us flrst consider the decoupling limit in which the expressions are simpler but the 
structure of the divergences is the same as in the full theory because higher powers (two or 
more) of A^^-^"^ generate flnite contributions only. In the decoupling limit Eq. f lA14p simplifles 

to 

- g (^6stml - 2 (5.'+3..rSA^^^rg;5,,, + 5..rg',^A^^^rg;,5,,+3,t) j ^ + finite . 

(16) 
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\ if 



9i+3 



9/ 



Qs 



a] 



b) 



FIG. 2: Decomposition of the squark-tadpole which is contained in diagram 2 of Fig. [T] as a 
subdiagram: In the decouphng hmit the squark tadpole is either proportional to one element /S. -- 
(diagram a)) or independent of A?-^^ (diagram b)). In the hrst case, it connects left-handed with 
right-handed squarks, while in the second case it is flavour and chirality conserving (proportional 
to 5 St)- The divergence of the piece proportional to Aj^^^ is absorbed by the counter-terms to 
and A'^ while the divergence of the piece stemming from diagram b) is canceled by a squark mass 
counter-term. 



Here we clearly see that in order to render the first term in Eq. ( fT6|) finite, which is flavour- 
diagonal (corresponding to diagram b) in Fig. |2]), a renormalization of the squark masses is 
necessary. On the other hand, for canceling the divergence of the second term in Eq. ( IT6|) 
(corresponding to diagram a) in Fig. [2]), which is proportional to A^^ , a counter-term to 
the Yukawa coupling and the A-teim contained in A^^^^ is necessary. The latter point can be 
seen as follows: In the decoupling limit the amputated chirality changing squark two-point 
function for q^, — )■ is given, at lowest order in a^, by 

-pii'-k \qRL-pjj' 

From this we can read off the common renormalization renormalization constant Zy of the 
Yukawa couplings the Ajj and the A^'^ terms, obtaining in the minimal subtraction 
scheme (DR or MS) 

Zy = 1-^2Cf-. (18) 

In fact, it turns out that this renormalization of the Yukawa couplings is necessary for 
maintaining supersymmetry with respect to the Yukawa coupling involved quark-quark- 
Higgs coupling and the one of the squark-squark-Higgs coupling. 
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C. Result for Cj/"^ retaining all powers of t^/MsusY 

For the Wilson coefficient C'j/"^ of the two-quark operator q^Pjiqi we write the general 
decomposition 

cjf = c-;,^"'" + ^ ^c<;' + (^) ' . («) 

From Eq. flA2p we directly obtain 

= |j (4 ^/.W^ts,.^^) • (20) 
Here we introduced the abbreviations 

Xt = m^Jmg, (21) 
and for later convenience we also define 

yt = mujmg,, Xf, = fi/mg, (22) 

where fi is the renormalization scale. 

According to the detailed description in the previous subsections, we decompose the 
Wilson coefficient Cj-/ into various pieces: 

^(2) ^ ^(2),1 ^ ^(2),2 ^ ^{2),3 ^ ^{2),4 _^ ^(2),^* _^ ^(2),shiftas _^ (-y(2),m^ ^™'5dr _^ (-7(2),T'P ^23) 

We freely made use of the unitarity of the mixing matrices and obtain 

C}?'' = E E (2 Wf, Wfl,, wit,,, Wf^3,s (2 -Ca)j^ ^ 

j=i s,t=i L U 



X 



X2) (1 -X?) 

(1 - x',r Li, (1 - xl) - (1 - x?)^ Li2 (1 - X?) + (x^ - x?)^ Li2 (1 - x?/xD 
—4 x^ (xj — x^) In(xs) In(xt) + 6 x^ {xf — 1) In(xs) — 6 x^ (x^ — 1) In(xt) 
+2 x2 (x2 - 1) In^xt) + 2 (xf + x^ - 3 x2 x2 + x^) ln2(x,)] } 

+ E (4 W^/i WAta, m-g Cf (2 - Ca) 
t=i I (1 - xt) 

X [(1 - 2x2) \n\xt) - 2 (1 - X?) In(xi)] } , 



(24) 
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x{A{l-xl){-xl + 2x\-l)U^{l-xl) 
-4(x^ + xtf {xs - xtf Li2(l - xl/xj) 
-4xf ((1 + xf) ln{xt) + 1 - x?) - 2 (1 - x2)2 ln(x,)) 
— ^ [48 (In(xt) In(xs) (x^ — — xj) + In^(xf) (x^ — x^) — In(xt) x^ — x^) 

+24 (In(xj) In(xs) x^x^ (1 + xf) + In(xs) x^x^ (x^ — 1) + In(xj) xf + In^(xj) x^) 

+12 (In(xf) x^Xj (1 — xf) + In(xs) x^ (1 — xf) — x'^xf) 

+6x2 + 6x2x^ + 30x^ + 18]}} , 

^(2),3 ^(2),2/ ^ N 

x{(-3CA + 2CF)Li2(l-x2) 

1 

^2{l-Xtf{l + Xtf 

X [trnj ((241n2(xt) - 12) x^ + (241n(xj) + 12) x?) 

+6 Ca ((3 \n\xt) - 11 In(xt) + 9) x^ + (3 In(xt) - 14) x? + 5) 
+3 Cf (-(2 In(xt) + 1) x^ - (12 \Yi\xt) - 12 In(xt) + 7) x? + 8 In(xt) + 8)] } } 
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f = E i -4 Wl, Wf^,, m, — ^ 

t=l \_ yL X 

X [tiuf {-2xi ((21n(xt) - l)x1 + 1)) 
+CA(3x?((21n(xt)-l)x? + l)) 

+ ^ (-(4 In(xt) - l)xj + (2 In(xt) + 3)a;2 - 8 In(xt) - 4) 



ln(x 



2\2 



(28) 



(2),shiftQ:s 
fi 



6 

E 



3fl 



E(ln(x,)+ln(2/,))-4C^ + 3C^ 



W^) 



(29) 



' ^MS "dr. 



CfF' = -2 m,Cl E { W^/. ^ts^^Y^ 

6 



xt 



2\n{xt)) (l-21n(x0+ln(x2)) 



j,j'=ls,t,s'=l L ^ ^ I'J-'-J 



2(2 In(x0-ln(x2)-l) a^^x^, In 



Xs^ 

X a 



+ In (xs) — x^, In [xgi] 



(x2 - x2,) (x2 - 1) (x2 - 1) 



In the MSSM we have 



3, C^ = 4/3, tr = 1/2 



and 



nt 



6, 



(32) 



To summarize: Eqs. f l20p and (1231) contain the full result for the Wilson coefficient Cjl"^ 
where the A-terms, the Yukawa coupling, the squark and the gluino masses of the MSSM 
are renormalized in the DR scheme, while gs stands for the strong coupling constant of the 
SM in the MS scheme, running with six flavours. The effective operators, or equivalently 
the Wilson Coefficients, are understood to be renormalized according to the MS scheme. 

So far, we discussed the derivations of Cj^"^. The corresponding result Cj^^^ for up- 
quarks can be obtained by replacing with and exchanging x and y. 
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D. Reduction of the matching scale dependence at NLO 



The purpose of our NLO calculation is also the reduction of the matching scale dependence 
of the effective Higgs couplings which can serve as an estimate of the theory uncertainty. 
This reduction is not only an improvement achieved by our NLO calculation, it also serves 
as an additional check of is correctness. 

As we will see in the next section, the quantity directly related to the Higgs couplings is 
Cj- defined as 

^jLR ^ ^jLRW + c%LRi2) ^ 1 (^cj^'cJ^RW + ^(D^f ^) + Q {a^ tan/3) . (33) 

We use in the following the decomposition Cjl"^ = Cjf^'^^^+C'jf^'^^^. At LO in our counting 

of as and tan/3 we have Cj^^'^^^ = Cj['^'^^\ 

C'fi^ (and thus also Cj^^^^^) at a fixed low scale /iio„ is obtained from Cjl"^ at the 
matching scale /iq via 

Cfj^^ (yUiow) = U (/iiow, Ato) C'f/^ (/io) • (34) 
This evolution is the same as for the quark masses in the SM. The explicit NLL expression 



can be taken e.g. from Eq. (4.81) in Ref. 32|. It is this expression which we use for the 



numerical study in Sec. IIV El when doing the evolution to the low scale fiiow 

However, for showing analytically the reduced matching scale dependence, it is sufficient 
to assume that the scale fiiow is close to the matching scale fio so that it is not necessary to 
resum large logarithms. In this case the evolution matrix U (/iiow,/Wo) can be expanded as 

(0) / 2 \ 



U ifiio^, /io) ~ 1 + «s (/io) In ) , i^^^ = QCf. (35) 



low 



At LO Cfj^^ depends only implicitly on the renormalization scale via the scale dependence 
of various parameters. For small changes of the original matching scale /iq to a new matching 
scale fi, we get 

C^r'^'if^) , ^^^^^\p, + S)\n(i4) . (36) 



The contribution involving /3o comes from expressing as(/i) in terms of ^^(yUo); while the 
one involving S is due to the corresponding manipulation of the squark- and gluino masses, 
the Yukawa couplings and the A (and y4')-terms. Together with Eq. ( I34l) and Eq. ( I35l) the 
variation of the matching scale leads to the following ratio 

U CJ'""' (f ) . 1 + fliM + s - 2^ 'lln I ^ I . (37) 




The explicit /i dependence proportional to in this ratio has to be compensated when 
going to NLO. 
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U(//iow,y")Qj W 

U(//iow,lTeV)C-j'^''(lTeV) 
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FIG. 3: Dependence on the matching scale /x of the one-loop and two-loop result for Cj^^ (fiioy^), 
using MsusY = 1 TeV and ^uiow = rn]y. Red (dashed): matching done at LO; blue (darkest): 
matching done at NLO matching. As expected, the matching scale dependence is significantly 
reduced. For the one-loop result, C'i^'^ is understood to be C'i^^^^^ (see text). 




The piece of C^j^^^"^^ with explicit scale dependence (with contributions from Eqs. ( 128|) . 
([29]), dS]) and from Eq. §^ through C^f and Cf^ ), can be compactly written as 



Air 



S -2tinf + 3Ca-^Cf + 



qLR{l) 



(/i) ln(x 



Using this information, we finally get at NLO 



U (/Xzot«,/io)C'f,-^^(yUo) 



1 + 



as expected. 



(0) 

/3o + 5-^-5 + 2trn/-3C^ + 3C^-^ 



2 In 



(38) 



(39) 



E. Numerics 



In this section we study the numerical importance of our two-loop corrections and the 
reduced matching scale dependence compared to the one- loop result. 

The matching scale dependence, as shown in Fig. |3] for SUSY masses of ITeV, is signifi- 
cantly reduced as expected from the previous subsection. Note that the relative importance 
of the NLO result is to a very good approximation independent of the size of A^- . 

The relative importance of the two-loop contribution to Cfj^^{fi) is shown in Fig. HJas 
a function of the matching scale For SUSY masses of 1 TeV the corrections lead to 
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FIG. 4: Relative importance of the two-loop corrections as a function of the matching scale /i. 
We see that the two-loop contribution is approximately +9% of the one-loop contribution for 
/i = MsusY = ITeV. C^f^ is defined in Eq. ([33]). 



a constructive contribution of approximately 9% compared to the one-loop result which is 
in agreement with Ref. 13|. Again, the relative importance of the NLO result is to a very 



good approximation independent of the size of A 



qLR 



F. Transition to the decoupling limit 

While the two-loop contributions calculated in this section are obtained in the approxi- 
mation p = rriq = 0, the results given in Subsection IIV CI still contain all powers v/M^vsy 
implicitly via the squark mixing matrices and the physical squark masses m^^ involved. 
The transition to the decoupling limit, in which all chirally enhanced corrections can be 
resummed analytically, can be done by the following prescription. 

In all parts of the genuine two- loop contributions listed above (Eq. f|24p -Eq. (1301] ) only 
two mixing matrices occur, except in Eq. and Eq. fl3ip . Eq. f l2^ contains the following 
combinations of mixing matrices and a loop-function / which depends on squarks masses 
rriq^ and rriq^ 

6 3 

E E ^t3,t wi:s,s w^ts.s fi^i ^t) ■ (40) 

s,t=l j=l 

Note that in the decoupling limit, the squark with index s in Eq. f HU]) must be a linear 
combination of right-handed squark only, since otherwise at least two chirality changes (two 
insertions of Af^-^^) would be necessary. Thus we can replace 

W,t3,, l^ts,. ^ = and xl ^ xl, , (41) 
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where k only runs from 1 to 3 and we defined 



2 

m 



L(R) 



The resulting expression 

E E Wf, W^,, Kif /(4„ x?) (43) 

can now be expanded in powers of f/MsusY which amounts at leading order to the replace- 
ment 

EU/'^U/'^* f( , \^ KdLL \dLR f^dRR fj--- i ^Im) ~ /(••• ? ^fln) 

t=l m,nj',j"=l 94 9n 

where the dots represent possible additional dependences on squark masses. Now we apply 
Eq. dS]) to Eq. (SS]) and use 

3 

E,d{LL)RR ,d{LL)RR. _ ,d{LL)RRr / . c-x 

■m fj ^^n ji — ^^m fi "rnn ■ {^O) 

i=i 

The final result for Eq. (HOj) in the decoupling limit is then 

E\dLL /\dLR\dRR f i^Rn^ ^Im) ~ f {^"ru^ ^Rn) fAr-\ 
^^m fj"^j"j' ^^n j'i ^2 _^2 " I'^^J 

m,n,j',j"=l 9m 9n 

For Eq. ( 13T]) a similar procedure works. It contains the following combination of mixing 
matrices with a loop-function depending on three different squark masses with the indices 
s, t and s' 

6 3 
s,t,s'=lj,j'=l 

(47) 

Note that the first term in Eq. (l47j) vanishes in the decoupling limit since it necessarily 
involves multiple chirality flips. For the second term two replacements analogous to Eq. ( 14T|) 
have to be performed and after using two times the relation in Eq. ( I45l) the decoupling limit 
of Eq. (HZD reads 



E\dLL \dLR\dRR f Lm^ ^ Lm^ ^ Rn) f Lim ^ Rm ^ Rn) /ao\ 
fj"^j"j' j'i ^2 _^2_ • y^^) 

m,n,j' ,j" = l <irn 'in 

This result involves the same combination of mixing matrices as the one in Eq. ( l46l) . 

To all other parts of Cji the rule in Eq. (jH]) can be applied directly in order to obtain 
the corresponding expression in the decoupling limit. 
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V. RELATIONS BETWEEN QUARK MASSES AND THE MSSM YUKAWA 
COUPLINGS AT NLO 



Beyond one-loop Eq. ffTTj) and Eq. f|T2|) for the determination of can easily be general- 
ized to higher loop-orders because the chirality changing self-energy (and also the resulting 
Wilson Coefficient) is still proportional to one element Afj^^ in the decoupling limit, as 
shown in Sec. IIVFI However, since we are dealing with order one corrections, we must spec- 
ify how we count contributions at higher loop-orders in a^. Cj/'^''^^ is proportional a^tan /3 

and Cjl'^^'^^ is proportional to a^tan/3. Here, tan/3 stands schematically for a chiral en- 
hancement factor, also including Aj^/ (Y^jMsusv)- We will count a^tan/? as order one and 
thus a^tan/3 as order a^- Since cj^^^^^ jg not chirally enhanced, the only relevant term 
in our approximation (of order «s) is the one-loop contribution. Thus, cj^^^^^ jg always 
understood to be the one- loop contribution proportional to a^. 

In order to derive the relation between the quark masses and the Yukawa couplings of 
the MSSM superpotential at NLO we also need to specify the renormalization scheme used 
for the matching procedure. Let us explicitly denote the renormalization scheme for the 
quantities in the matching condition Eq. flTT]) (at the scale t^susy) which is important at 
NLO: 

v,Y-^' = mf= (.,Y-r + q'^''' + Cl'^''') X (l + i (Cf' + Cr^)) . (49) 

Again, Y^^l^^ is the Wilson Coefficient induced via the Yukawa coupling of the MSSM. This 
means at the matching scale it is given by: 

Y,lT' (/^susy) = (/isusY) = (l + ^C^) (/^susy) • (50) 

In our counting in ag and tan /3 the renormalization scheme for Cf^^^ and Cj^^^ is irrelevant. 
Note that the quark mass m^^ is understood to be evaluated at the matching scale. Further, 
one should recall from the last section that despite the fact that we renormalized Cf^^^ in 
the MS-scheme, it contains parameters given in the DR-scheme, e.g. Y^^ = Y^. Since we 
are jnterested in Y^\ the Yukawa coupling of the MSSM superpotential, we must express 
y^fee i'^ Eq. (H9|) in terms of Y^ via Eq. ( |50|) so that we can solve for Y^. 

In conclusion we arrive at the NLO generalization (order a^tan/3) of Eq. fITT]) : 

™MS /ydLR 

- „,(i + ^c, + ta. ' ' 

with C^j^^ defined in Eq. ( I33l) and the corresponding equation for ef. Here e^^^"* and ef''^-' 
are defined in direct analogy to Eq. (1131) . Further, the Wilson Coefficients appearing here 
are assumed to be in the decoupling limit. Eq. ( ISTl) constitutes the NLO determination of 
the Yukawa coupling of the superpotential. When inserting later the Yukawa coupling into 
the Wilson Coefficients, one has to use this relation^. 



The generalization to the CKM matrix can be achieved following the procedure of 0, Q 
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The electroweak contributions (involving charginos and neutralinos) to the relation be- 
tween the quark masses and the Yukawa couplings are in most regions of the parameter space 
subleading compared to the strong contributions. However, the LO electroweak corrections 
are easily as large as the NLO SQCD corrections and should be included in a numerical 
analysis. This can be achieved by simply adding the corresponding contributions to C'^^^ 
and if in Eq. (IHTj) . 



VI. EFFECTIVE HIGGS VERTICES 



In order to derive the effective Higgs-quark-quark couplings^'' we have to assume that 
the external momenta (flowing through the Higgs-quark-quark vertex) are much smaller 
than the masses of the virtual SUSY particles running in the loop. This assumption limits 
the applicability of the resulting Feynman rules. If m^^o , m^o , m//± -C Msusy {H^,A^,H^ 
denote the neutral CP-even, CP-odd and the charged Higgs boson, respectively), the ef- 
fective Feynman rules can be used for the calculation of all fiavour-observables (also if the 
Higgs is propagating in a loop) and for processes with a Higgs on the mass shell. If the 
hierarchy mj^^o , m^o , m/^± ^ Msusy is not satisfied the effective Higgs vertices can still be 
used for processes in which the momentum-flow through the Higgs-quark-quark vertex is 
small compared to MgusY which is true for all low-energy flavour observables with tree-level 
Higgs exchange (like Bd^s /^^/^~) t^v oi the double Higgs penguin contributing to 

AF = 2 processes). 

As discussed in the introduction we use an effective field theory approach in our study of 
the Higgs-quark-quark couplings which simplifies the calculations significantly. This means 
that we match the MSSM on the 2HDM of type HI at the scale Msusy rather than calculating 
the Higgs-quark-quark coupling within the MSSM. 

Let as first consider the effective Lagrangian of a general 2HDM (including Higgs-quark- 
quark couplings and kinetic terms): 



+ df Ri0 [Sfi - Rjr) diR + UfRi0 {5f, - R]r) Ui R 



In principle also the renormalization of the Higgs potential should be addressed. Our derivation of chirally 
enhanced flavour effects does not depend on the specific relations between Higgs self-couplings and their 
masses. Since no chirally-enhanced effects occur in the Higgs sector, it is consistent to use the tree-level 
values for the Higgs parameters. However, one can as well use the NLO values for the Higgs masses and 
mixing angles which might be even better from the numerical point of view. 
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where adding the hermitean conjugate of the terms involving Higgs fields is implicitly meant. 
The Higgs doublets are defined as 



(53) 



In Eq. (1521] a, b denote 5'f/(2)j^- indices and is the two-dimensional antisymmetric tensor 
with ei2 = —1. We introduced the holomorphic couplings Ej^^, the non-holomorphic cou- 
plings EJ^"" [q = u, d) and the contributions to the kinetic terms R^f^^ and Here the 
superscript "ew" refers to the fact that these terms are given in a weak-interaction eigenba- 
sis. In Eq. ( I52l) we already anticipated the MSSM where the terms E^^]^'^'^^ , and I^^^^ 
are loop-induced but Yf^l^ and Yf^l^ are generated at tree-level via the MSSM Yukawa 
couplings^ ^. 

In order to connect the effective theory to the MSSM we go to the super-CKM basis, in 
which the Yukawa couplings are diagonal, by rotating the fields 

^,L,K^U]l^^''^'\^L,R. (54) 

such that 

Tj-q L{0]*^qtmej-rq Bio) _ ^q, r /rrN 

^kf ^kjcw^ji — -' trcc'^/j ■ {'->'->) 

We now break the electroweak symmetry and write the effective Lagrangian in component 
form: 

+ df ^V^f^ ( [Y^iJ,, + El) Hl^ - E'-Hj) n 
-dfj (Y^Jf^ + Ef^ + E%hA d, n 



-UfL{ {YZSf^ + El) Hl^ + E'lH}) n^ R 

+ df Ri0 {6 fi - Rji) diR + UfRiB [5fi - R^) Ui r 
+ df Li0 {Sfi - Lj-) diL + U f LiB {5fi - LJ-) ui l 
-djL[ [Y^^iJfi + Ej^i V, + Ef^^) d. R 
-UfL{ [Y^JJj, + El) + E'lvd) Ui R , 



(56) 
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In principle, without knowing anything about the MSSM, the holomorphic corrections could be absorbed 



into an effective Yukawa coupling (and also the corrections to the kinetic terms Rj^^ ™ and would not 
be physical). However, once we go back to the MSSM with the SUSY breaking terms as input parameters, 
also the holomorphic corrections become physical. 
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where V^^^i = f/« i(o)tf/dL (o) jg physical CKM matrix, but rather the CKM matrix 

generated by the misahgnment of the Yukawa couphngs. Adding the hermitean conjugate 
of the mass terms and the terms involving Higgs fields is tacitly understood. The terms 



fi kf kj ji 

j^iq _ jjQ L{0)* j^/q cw jjQ R{0) 
fi kf kj ji 

Rj^ = Ul^^^'^^RlYU"/^''^ (57) 

Td _ Tjd L{Qi)* J qcwjjd L{0) 
^fi - ^kf ^kj ^ji 

jju _ jju L{0)* j^q cw jju L{0) _ j^d y W* 

fi kf kj ji fk kj ij 

are now given in the super-CKM basis. Note that this is the same basis as the one in which 
the effective Lagrangian of Eq. ([2]) is given (and the same basis in which we calculated the 
MSSM contributions to the Wilson Coefficients). Thus, comparing the last four lines of 
Eq. (1561) to Eq. ([2]) we have the following relation between the Wilson Coefficients and the 
terms of the 2HDM III Lagrangian (at an arbitrary loop order): 

^d ^ '-'fiA ^,d _ '-'fi 



Vd ' ' 

Er = ^, E'fl = ^^^, (58) 

^ Vy, ^ Vd 

rq ^ryqLL jyq y^qRR 

^fi - '-'fi ) ^fi - '^fi ■ 



Now we want to go to the physical basis with flavour-diagonal mass terms and canonical 
kinetic terms. As a first step we render the kinetic terms canonical by a field redefinition: 

1 



qn^ [^ij + -L%j Qj L , 
QiR^ i^ij + \Rtj] QjR- ^^^^ 



Consider now the quark mass matrices. The redefinition of the fields in Eq. (159!) also leads 
to a shift in down-quark mass matrix so that it is now given by 



fi ~ \ fi + fi J Vd^fi - '-'fi + fi ■ 
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where we have defined 



E 



fi 

^ fi 
q LR 



1 ^ 

2 



1 3 



rqtrce _ r , / r^qLLyrqi r , -t^g/ r ^qRR 



(61) 



C 



fi 



c 



q LR 
fi 



El ^q LL^q LR . ^q i^n^q r 
^ . . \^fj ^ji + ^fi ^ji 



yq LR^qRR 



Note that the quantities with a double hat contain also the contributions from fiavour- 
changing LL and RR Wilson Coefficients, while the quantities with one hat (see Eq. ( 1331 
and Eq. ( 1671) ) only contain the ffavour-conserving LL and RR Wilson Coefficients. 

We now diagonalize the quark mass matrices by a bi-unitary transformation 



'kf 



q."fi 



(62) 



where the rotation matrices 
/ 



-m 



13 



1 


q 

mh 


m?3 


\ 










1 


q 

"^23 








nig^ 




q* q* 


q* 

-"^23 


1 




777,^2 TTT'q^ 




/ 



R 



m 



q* 
21 



m 



q* 
31 



m. 



12 



m, 



-m. 



21 



m 



93 
32 



-m 



31 



+ 



92 

minm 



93 



32"^21 



32 



93 



(63) 

are obtained from a perturbative diagonalization of the quark mass matrix^^. 

Switching to the physical basis in which the quark mass matrices are diagonal, these 
rotations modify the effective Lagrangian as follows 17|: 



reff _ ^ Tju L*T/(0) 



^Hf - E'^,^ {Hi + tan (/3) Hf) \ R 

Vd 



-J- rl Tr ^*V^ 



k'k 



- df LUt/-' I '-^m* + Eg (/72 - tan (/3) H'/) ] Uf.^d, « 



^fL^kf 




(64) 



'^^^ Hl^ + (i/] - cot (/3) i/f ) ) f/" ^u, n + h.c. 



^kj 



where we skipped the mass terms and the kinetic terms. This can be further simplified by 
using that the physical CKM matrix is given by 



fi 



(65) 



^"^ Note that these rotations are identical to the ones obtained in the diagrammatic approach (see Ref. [17 1 
for details). 
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In addition we define the abbreviations 



( 



-^22 '-12 



fq LR j^/q ( LR _ LRpff LR 



33 I '-13 





Piq ^q LR piq LR 

-^22^21 -^^33^23 

Piq ( -qLR _ p.q LRp_q LR\ piq -g LR 

-^33 I "-31 '-32 '-31 j -^33'-32 ^ 



12 "-23 



(66) 



fi 



Note that in this expression only quantities with a single hat defined as 



E 



E 



LffLljf- + Lljf- K 



and cl^^ defined by combining Eq. fl33|) with 



{')|J 



fi Si 



(67) 



^g LR 



gLi? 



max{mg^.,mgj 



(68) 



enter. This is in agreement with the finding of Ref. 20| that the effect of the flavour- changing 
LL and RR self-energies drops out in the effective Higgs vertices. 

Finally, in order to arrive at the effective Feynman rules we project the flelds and 
onto the physical components , and as 

Hi = ^ (iJ° sin a + h^cosa + iA^ cos /?) , 
v2 



V2 

Hi* = cos{P)H-., 
Hi = sin {/3)H~. 



H^ cos a — sin a + iA^ sin /3) 



(69) 



Using Eq. (!65|) . Eq. (!66|) and Eq. ( 169|) . the effective Lagrangian in Eq. (IMIl leads to the 
following effective Higgs-quark-quark Feynman rules^^ shown in Fig. [5] (note that the CKM 



Note that some of the Higgs-quark-quark couphngs arc suppressed by a factor cos /3 or sin a stemming 
from the Higgs mixing matrices. If one decides to keep these suppressed couphngs, one should be aware 
of the fact that they receive proper vertex corrections in which the suppression factor does not occur and 
which arc thus tan /3-cnhanccd with respect to the tree-level couplings. Such enhanced corrections to the 



coupling of H to right-handed up-quarks are important for 6 — ^ 57 33|, |3J| . 
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I 
I 



H- 



(If 



di 



Uf 



* U Ufd, -TR + i Ufdi 



FIG. 5: Higgs-quark vertices with the corresponding Feynman-rules. The couphngs with exchanged 
chirality structure are obtained from Eq. ([70|) by using ^^^q^ = ^q^^*- 



matrix V in the charged Higgs coupling is the physical one): 

Tufu? = xt f ^5/. - E'- cot (3) + x'/E';, , 



J2 cos /3 V;* (^^5,-, - E';, tan , (70) 



where for = {H^, /i", A") the coefficients are given by 

Xj = ^cosa, ^^sma, ^^smp , x , = ^sma, ^cosa, ^^cosp 

\ V2 V2 V2 J \ V2 V2 V2 

(71) 

It is important to keep in mind that the c^j in Eq. ( 166|) must be calculated using the 
quantities and V^'^^ of the MSSM superpotential. 

Note that without the non-holomorphic corrections E'^ the rotation matrices W^^'^ would 
simultaneously diagonalize the effective mass terms and the neutral Higgs couplings in 
Eq. (^^. However, in the presence of non-holomorphic corrections this is no longer the 
case and apart from a flavour- changing non-holomorphic correction also a term proportional 
to a flavour- conserving non-holomorphic correction times a flavour- changing self-energy is 
generated. 
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A. Effective Higgs-quark-quark vertices at NLO 



The effective Higgs-quarlc-quarlc vertices at NLO in the MSSM are obtain in the following 
way: After inserting the definition for E'j- (see Eq. (l66l) ) into Eq. flTOj) we express E^^"^ 

through C^rj'^ according to Eq. (158|) . 



VII. CONCLUSIONS 



In this article we computed the genuine two-loop SQCD corrections to the chirality- 
changing quark self-energies. In the limit where the external momentum and the quark mass 
are zero, we presented relatively simple analytic results without making further assumptions 
on the SUSY spectrum. Due to the one-to-one correspondence (in the decoupling limit) 
between chirality-changing quark self-energies and Higgs-quark-quark vertices, this is an 
efficient and elegant way of calculating at the same time not only effective Higgs vertices, 
but also the Yukawa couplings and CKM elements of the MSSM superpotential in terms of 
the physical quark masses and the physical CKM matrix. 

Our next-to-leading order results increase the values of Wilson Coefficients C of the 
operators qf Pr Qi by approximately 9% compared to the values obtained at leading order . 
At the same time the matching scale uncertainty of the effective Higgs-quark-quark couplings 
and of the corresponding Wilson Coefficients is significantly reduced (see Fig. [3]). 

We resummed all chirally enhanced corrections modifying the relation between the quark 
masses and the Yukawa couplings of the MSSM superpotential up to order a""'"Han"'/3 (see 
Eq. ( 151]) ). The resulting MSSM Yukawa couplings can be used for a precision study of 
Yukawa unification. Furthermore, using these Yukawa couplings, we derived effective Higgs- 
quark-quark vertices (see Eq. (fTOj) ) entering the calculation of FCNC processes and also of 
Higgs decays, as long as the momentum transfer is small compared to the SUSY scale. 
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Appendix A: One-loop results 



Here we summarize various one-loop results necessary for the two-loop calculation of the 
chirality flipping self-energy (see 35| for details). Unless stated otherwise, all expressions 
appearing in this appendix were obtained in dimensional regularization. The matrices 
diagonalize the squark mass matrices according to Eq. ([6]) and we use the definitions: 



m-g 



Trig 



X,. 



m-g 



(Al) 
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1. Self-energies 



Here we give the explicit one-loop results for quark, gluino and squark self-energies in 
dimensional regularization, where we put D = 4 — 2e and write the renormalization scale 
in the form fie'^^'^ /{\/Att). Our conventions are such that the calculation of the truncated 
self-energy diagrams give — zS. 



a. Quark 



The one-loop quark self-energies induced by gluinos and squarks are given by 



x2 (In^ {xD - 2 In {xD - 2 l n(a;^)ln (a:^)) ^ ^ 



2 1 - x2 



-^WlW!:C,(^l + ln{xl) + 



2(l-a:2) 



+ 0(e) . 



(A2) 



Using unitarity, we can replace i?o(0;m?,m?J by [-Bo(0; m?, m?J — i?o(0;m?,0)] in the first 
line of S^^^^(O). This we did when writing the explicit expression. 
The ordinary gluon correction reads in Feynman gauge 

j^,LL,RB. (p2) ^ _ 2) B, [p^- ml 0) , 

= —Cf — + ^ In — - ^ + In ^ - 1 6h 



+ 0(£) 



S?/,f '^^ (p2) = ^CFd m,^ Bo {p'; ml , O) 5^, 



(A3) 



Up 

TT 

+ 0(6 
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Note that Eq. (1A3P and Eq. (lA4p are given in dimensional regularization (not dimensional 
reduction). 



b. Gluino 

Here we assume that of the three gaugino masses the gluino mass is chosen to be real 
which is always possible. For the gluino self-energy the part induced by a gluon reads 

^9 (P') = ^C'a {dnigBo {p^; m^, o)+p{d-2) B, ^1 O)) , (A5) 

which decomposes for on-shell gluinos into 



K)=-^C^Q + 2 + ln(x3)+0(.), 



(A6) 



■^g LL,RR (^2 
^99 



where we inserted the explicit expressions for the loop-functions. The part of the gluino 
self-energy with squarks and quarks as virtual particles in the approximation ntq = is 
given by 

j,,|i.H,RL = 0, 

Sir™!?^) = ^2(rf (B,(/i0,m?)+ft(p^0,O), (A7) 
where the latter reads explicitly for on-shell gluinos 



^ + 2 + In (xy ) 



+ 0(6) , 

(A8) 

with Tij = 6. The quantities Il^^{m^) and S|'^(m?) which appear in eq. (1A26P are defined 

as 

Sf '^«(m|) = S?/^'^« {ml) + Eg ^^'^^ (m|) . (A9) 



c. S quark 

For the squark self-energy we have 
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where the parts refer to the squark self-energy with gluon 

a. 



An 



Cf (2 + ml) Bo {p'; ml, O) - Ao {ml)) 6st , 



K) = 3 ^C,ml (i - In (.1) + In (xj) + 5., + O 



the squark self-energy with quark and gluino 



iP') = (Ao {m^ + {ml - Bo (p^ m?, O)) 5., , 



a.. 



QsQt \ Qsl 271" 



+ 3-2x2 + (2-x2)ln (xy 



1 



- 2 In (1 - - iO) 



5.* + O ie) 



(All) 
(A12) 



(A13) 



and the squark tadpole self-energy of Fig. [2] (for up (down) type squarks only the diagram 
with internal up (down) squarks is non-zero): 



yqq _ 



^C, {6.,Ao {ml) 



3 6 



-2 E E [w^:3sWUss'w;:w^,+w!:wi,w^i,,wj^^^^ 

i,j=ls'=l ^ ^ ^ 



^ + l-ln(x2)+ln (x2)^ 



3 6 



i,j=ls'=l ^ 



X I - + 1 - In (x^,) + In (x^) ) +0 (e) 



(A14) 



Note that S™-^ is independent of the external momentum. The part proportional to Sst in 



eq. (1A14P is due to diagram b) of Fig. [2] while the second part, which is proportional to at 
least one element A^^-^"^, is generated by diagram a). 

Note that in the sum of all contributions to the diagonal squark self-energy there is no 
divergence proportional to and thus no wave-function renormalization is needed in order 
to render the diagonal squark two point function finite. 



2. Loop-functions 

The one-loop functions Ao{m'^), Bo{p'^]m\,m^ and i?i(p^;m^,m|) in the previous para- 
graph are defined as 

2, _ 167r2 /i^^e^^ f d'^i 1 



Aoim' 



i {AttY J {27rY P - m2] 



(A15) 
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Bi (p^; ml, ml) 



i {Any J {2txY - ml] [(£ + pf - ml] 



(A16) 



(A17) 



i {4ny J {27iY - ml] [(£ + pf - ml] 

The function BQ{ml,ml) which also appears, is an abbreviation for Bo{0]ml,ml). We give 
now relations among these functions and explicit versions for specific arguments 



Aq [m?) = m^ 
Bo {ml, ml) = 



m^ 



12 



m^ 



2 _ ^0 {ml) - Ao {ml) 



Co {ml,ml,ml) 



1 1 
mf — m2 



dBo {ml, ml) 



1 



5o(p^;m^O) = --ln — 



m 



1^' 



1 2 

m — p 

p2 



In 



m? — — iO 



m^ 



Bi (p2;m2,0) 
Bi (p2;0,m2) 



Bi{0; ml, ml) 



with X = m2/mi. 



2p2 
1 

2^ 



[Ao{m'^) - (p2 + m^) Bo {p^; m^, 0)] , 
2 [-Ao{m^) - (p2 - m^) Bo {p^; m^ , 0)] 



0{e), 



1/1 , . 3 - 4a;2 + + (4x2 - 2a;^) In (x^; 

'2 + — 



2(1 



(A18) 



3. One-loop renormalization and counter-terms 

a. One-loop counter-term diagrams 

Squark-mass counter-term diagram: 

- (y ^ - - 

^f.r"^ = ^Cpm-, E 5mlW%Wt;,sCo ( 



=1 



2^ 

'-C^m.EWf.Wf*,,^ 



2 2 



27r 



(1-x: 



[In (x2) + l-xl 



..(ln(.S(l-.^.ln(.D)4ln^(.D.(ln(.D-l).^.l 



(A19) 
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Gluino mass counter-term diagram: 



a 
2n 

a 

2vr - s=i 
—e 



6 



= T^C'p^™, E W},Wf:,, {Bo (m?,m|) + 2m|Co ( 



)) 



2 1-x 



tts counter-term diagram 



2^2 



[4 (1 - x^) + 2 (1 + ln(x2) - (1 + xl) - xl) 
+ {A{l-xl) + 2{l + xl) Hxl))Hxl)]] 



yL/i as CT — f 

(5a 



s=l 



1 

i+3,s 



x1 In {x^ 



2(l-x2) 



(In^ (a;2)-21n(x2)-21n(x2)ln (x^)) 



(A20) 



(A21) 



b. Renormalization of the Yukawa couplings in the MSSM 



Due to supersymmetry, the renormalization of the Yukawa couphng in the quark-quark- 
Higgs vertex and the one in squark-squark-Higgs vertex Y^^ must be identical^^. Indeed, 
we exphcitly find that the counter-terms for these couphngs are the same 

ygj,9i{0) _ _|_ ^ygi.gi ^yqi.Qi _ _^'^(J^Y'^"^^ (A22) 

An e 

which even holds in the MS— scheme and in the DR— scheme at the one-loop level. 

c. A-term renormalization 

In the approximation rriq = the SQCD renormalization of the A-terms is the same as 
of the Yukawa coupling^^. 

d. Squark mass renormalization 

We write the connection between the squares of bare and the renormalized squark masses 

as 

K)' = K-J' + K.- (A23) 



This also includes that the renormalization of the Yukawa coupling entering the squark mass matrices is 
the same as the renormalization of the quark-quark-Higgs coupling. 

If TOg 7^ the quark-gluino correction to A-terms induced flavour-non-diagonal (divergent) corrections. 
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From Eq. (1A12I) . Eq. (1A13P and Eq. (1A14P and by taking into account that the second term 



of Eq. (]A14p only renormahzes the Yukawa couphng (and the A, A' terms), we can easily 



read of Sm'^^. We obtain in the MS— scheme: 



6ml = ^C,ml ((x? + 4)-x?) 1, (A24) 

where the contribution proportional to {xf + 4) comes from Eq. (1A12P and Eq. (lAlSp while 
the term —xf stems from the part of Eq. f lA14p proportional to 6st- 



e. Gluino-mass renormalization 

We decompose the gluino self-energy according to Eq. ([1]). Expressing the bare mass 
(marked with the superscript (0)) in terms of physical one 

= nig + Srrig , (A25) 

we get in the on-shell scheme 

6mg = -m-g Sf (m|) - Sf (m|) . (A26) 

For details see Ref. [35j]. In the MS-scheme only the divergence of the right-hand side enters, 
i.e. we get in this scheme 

a 1 
Smg = -^rrig {3CA-2trnf) -. (A27) 

/. Renormalization of gs in the MSSM 

In lowest order, the strong coupling constant involved in Cjl"^ is Yukawa type. The 
relation between the bare and the renormalized version reads ^f^y = (1 -|- 5Zg^Y)9s,Y^ where 
the renormalization constant in the MS-scheme is given by 

(A28) 

Note that at one-loop the renormalization constant is the same for the MS-scheme and the 
DR-scheme. 



=4^ 



truf 



5c, 
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